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Abstract. We present a detailed gravitational mass measurement based on the XMM-Newton imaging spectroscopy 
analysis of the lensing cluster of galaxies CL0024+17 at z = 0.395. The emission appears approximately symmetric. 
However, on the scale of r ~ 3.3' some indication of elongation is visible in the northwest-southeast (NW-SE) 
direction from the hardness ratio map (HRM). Within 3', we measure a global gas temperature of 3.52 ±0.17 keV, 
metallicity of 0.22 ± 0.07, and bolometric luminosity of 2.9 ± 0.1 x 10 44 /if 2 erg s _1 . We derive a temperature 
distribution with an isothermal temperature of 3.9 keV to a radius of 1.5' and a temperature gradient in the 
outskirts (1.3' < r < 3'). Under the assumption of hydrostatic equilibrium, we measure gravitational mass and 
gas mass fraction to be M 20 o = 2.0 ± 0.3 x 10 14 /i^Mg and / gas = 0.20 ± 0.03h~g /2 at r 20 o = 1.05h~g Mpc 
using the observed temperature profile. The complex structure in the core region is the key to explaining the 
discrepancy in gravitational mass determined from XMM-Newton X-ray observations and HST optical lensing 
measurements. 
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1. Introduction 

One of the optically most prominent, but also most puz- 
zling distant lensing galaxy clusters, is CL0024+17 at a 
redshift of z = 0.395 (Gunn & Oke 1975). Its high galaxy 
density led to the early discovery by Humason & Sandage 
(1957), and the cluster has since been the subject of many 
studies. It was one of the first clusters to display the so- 
called Butcher-Oemler effect (Butcher & Oemler 1978; 
Dressier & Gunn 1982; Dressier et al. 1985), and Schneider 
et al. (1986) described it as a very rich optical cluster. 

Koo (1988) discovered gravitational arcs in 
CL0024+17, and it was subsequently studied exten- 
sively (Mellier et al. 1991; Kassiola et al. 1992; Kassiola 
et al. 1994; Wallington et al. 1995; Colley et al. 1996; 
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Smail et al. 1997; Tyson et al. 1998; Broadhurst et al. 
2000; Shapiro & Iliev 2000; Treu et al. 2003; Kneib et al. 
2003). A total of eight arc-like lensed images from the 
same background galaxy were identified (e.g. Tyson et al. 
1998). The redshift of the lensed galaxy was determined 
by Broadhurst et al. (2000) to be z = 1.675. Bonnet et 
al. (1994) detected a weak shear signal in this cluster 
out to a radius of 2.1 h^Q Mpc and inferred a high grav- 
itational mass of about 1.4-2.9xl0 15 fc^Mg. Recently 
Kneib et al. (2003) revised this lensing mass estimate 
downward on the basis of a very detailed Hubble Space 
Telescope (HST) mapping of the cluster which yielded a 
gravitational mass of M 2 oo = 5-f-i.o x lO 14 /if o 1 M out 
to r2oo = l-7h^ Mpc. In spite of the massive appearance 
of CL0024+17 at optical wavelengths, the cluster is 
relatively faint in X-rays with a ROSAT determined 
luminosity of L x = 1.22 ± 0.08 x 10 44 ft.f 2 erg s" 1 in 
the 0.1-2.4 keV band (Bohringer et al. 2000), an ASCA 
determined temperature of 5.71" 4 ) \ keV (Soucail et al. 
2000), and a Chandra temperature of 4. 47^^27 ke ^ 
(Ota et al. 2004; M 200 = 4.6±^ x lO 14 /if o 1 M to 
?"2oo = Mpc). In addition to the surprisingly low 

X-ray luminosity in comparison to its optical prominence, 
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Fig. 2. Chandra image in the 0.5-2.5 keV band (left) and surface brightness profiles (right) in annuli using the 0.5-2.5 keV 
image with azimuths 30-120° (counter clockwise) from north to east (black) and from south to west (grey), respectively. 



the detailed analysis of the X-ray observations yield a to- 
tal cluster mass in the ranee M 200 = 2-4.6 x 10 14 h^M Q 
that is 1.3-3 times lower than the weak lensing deter- 
mined mass (Soucail et al. 2000; Bohringer et al. 2000; 
Ota et al. 2004). The discrepancy is even larger at the 
arc radius ~ O.USh^ Mpc (e.g. Broadhurst et al. 2000). 

A hint of the explanation for this mass discrepancy 
came from the detailed analysis of the galaxy dynamics 
in CL0024+17 based on ~ 300 redshifts of galaxy mem- 
bers (Czoske et al. 2001; 2002). Czoske et al. (2001; 2002) 
found that the line-of-sight velocity distribution is not that 
of a relaxed cluster and is at least bimodal. They also 
demonstrated that the redshifts can approximately be ex- 
plained by a line-of-sight merger of two systems with a 
mass ratio of the order of 1:2. Also the weak lensing anal- 
ysis shows a mass distribution with substructure, modeled 
as a bimodal distribution of two systems with slightly dif- 
ferent central positions in the plane of the sky (Kneib et 
al. 2003). The infrared observations imply significant star 
forming activity with Star Formation Rates (SFRs) one 
to two orders of magnitude higher than those computed 
from the optical. The underestimation of the SFRs in the 
optical is due to absorption by dust (Coia et al. 2003). 

Using the capability of XMM-Newton to perform 
imaging spectroscopy at high angular resolution and thus 
allowing us to study the density and temperature struc- 
ture of the intracluster medium (ICM), we performed an 
XMM-Newton observation of CL0024+17 to shed new 
light on this enigmatic system. 

This paper is structured as follows. In Sect.|2 we de- 
scribe a double background subtraction method, which is 
developed to provide precise spectral background removal. 
In Sect. we analyse the properties of the hot gas in the 



galaxy cluster CL0024+17, then determine the total mass, 
projected gravitational mass, and gas mass fraction based 
on precise temperature and gas density profiles. In Sect. 01 
we discuss the structures in the cluster and possible solu- 
tions for the discrepancy between the X-ray and HST lens- 
ing measurements in the gravitational mass. In Sect.0 we 
present our conclusions. We adopt a flat ACDM cosmology 
with density parameter f2 m = 0.3 and Hubble constant 
H = 70 km s- 1 Mpc" 1 . Thus 1' = 0.320/if x Mpc. All 
coordinates are given in epoch J2000. We adopt the solar 
abundance values of Anders & Grevesse (1989). Error bars 
correspond to the 68% confidence level, unless explicitly 
stated otherwise. 



2. Method 

2.1. Data preparation 

CL0024+17 was observed on Jan 6th, 2001 for a to- 
tal exposure of 52.1 ks, 52.1 ks and 48.3 ks for MOS1, 
MOS2 and pn, respectively, by XMM-Newton with the 
European Photon Imaging Camera (EPIC) in standard 
Full Frame (FF) mode and Extended Full Frame (EFF) 
mode for MOS and pn, respectively. For all detectors, 
the thin filter was used. Data reduction and calibration 
was carried out using the XMM-Newton Science Analysis 
System (SAS5.4.1). The central position of the observa- 
tion is R.A. = 00 ft 26 m 35 s .7, decl. = 17°09'35".8. 

Above 10 keV, there is little X-ray emission from the 
cluster due to the low telescope efficiency. The particle 
background therefore completely dominates. The cluster 
emission is not variable, so any spectral range can be used 
for temporal variability studies of the background. The 
10-12 keV (12-14 keV) energy band (binned in 100 s in- 
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Fig. 3. HST mosaic image of the X-ray center of CL0024+17. Superposed X-ray contours of XMM-Newton as shown in Fig. 
indicate the elongation in the NW-SE direction on the scale of r ~ 3.3' (X-point). 



tervals) was used to monitor the particle background and 
to excise periods of high particle flux for MOS (pn). In 
this screening process we use the settings FLAG = and 
PATTERN < 12 (PATTERN < 4) for MOS (pn). 



We reject those time intervals affected by flares in 
which the detector countrate (ctr) exceeds a threshold 
of 2er above the average ctr, where the average and the 
variance have been iteratively determined from the ctr 
histogram below the rejection threshold. We screened the 
data using the thresholds of 21.4, 22.3 and 56.2 for MOS1, 
MOS2 and pn. The net exposure time is 46.6 ks, 46.0 ks 
and 42.9 ks for MOS1, MOS2 and pn, respectively. 



2.2. Background analysis 

We do not expect any cluster emission outside 5.4', the 
virial radius estimated from the redshift (z=0.395; Gunn 
& Oke 1975) and the Chandra average temperature mea- 
surement (4.47 keV, Ota et al. 2004). Thus the cluster 
emission covers less than half of the field of view (FOV) 
of the XMM-Newton telescope detectors. Actually the sig- 
nificant cluster emission extends only to 3' in the surface 
brightness of the XMM-Newton observations. Thus we use 
the 5 < r < 5.5' region in detector coordinates to compare 
the background conditions of the XMM-Newton pointing 
of the Chandra Deep Field South (CDFS, which we use 
as background field for the spectral analysis) with that 
of CL0024+17. The background of CL0024+17 is higher 
than that of the background field (CDFS) by - 30% for 
MOS and ~ 10% for pn on average in the complete en- 




Fig. 6. Projected temperature (left) and metallicity (right) profiles measured from the XMM-Newton data in the 1-7.8 keV 
band (solid, black) and 0.4-7.8 keV band (solid, grey) with point sources excluded. A solid curve in black indicates the best 
Gaussian fit of the temperature measurements in the 1-7.8 keV band with confidence intervals as dotted curves. Additionally 
the similarity temperature profiles in Markevitch et al. (1998; grey) and in Zhang et al. (2004b; black) are indicated by solid 
outlines. Our Chandra temperature measurements of CL0024+17 (dashed, grey) and the temperature profile of Sersic 159—03 
(grey shadow; 1 a interval) in Kaastra et al. (2001) scaled to CL0024+17 are presented for comparison. 




Fig. 1. Flat fielded, point sources masked and adaptively- 
smoothed XMM-Newton image of CL0024+17 in the 0.5-2 keV 
band with logarithmically scaled contours. Additionally we 
show the galaxy pair (cross) in Czoske et al. (2001; 2002) and 
the substructure (X-point) in kneib et al. (2003). 



ergy band, in which the difference is mainly produced in 
the low energy band below 1 keV. 

In the spectral analysis, the regions not affected by 
cluster emission, enable us to study the residual back- 
ground in the source observation compared to the back- 
ground observation. We applied a double background sub- 
traction method as described in Zhang et al. (2004a) using 
the results of the residual background modeling. In sum- 
mary the spectral analysis is performed in two steps using 
the software package XSPEC11.3.0: (i) A power law model 
for the residual background (background difference) is ob- 
tained in XSPEC from a comparison of the outer region 
(5-5.5') of the target and background fields (cf. Table 
(ii) The spectral modeling is performed in XSPEC with 
the cluster region of the target field as source data, the re- 
gion in the same detector coordinates in the background 
field as background and the residual background as a sec- 
ond model component with model parameters fixed to the 
values found in step (i). The overall spectra were fitted 
by a "mekal" model (Mewe et al. 1985; Mewe et al. 1986; 
Kaastra 1992; Liedahl et al. 1995; Arnaud & Rothcnflug 
1985; Arnaud & Raymond 1992) with Galactic absorption 
(Dickey & Lockman 1990; n H = 4.23 x 10 20 cm" 2 ). The re- 
sults with and without introducing a residual background 
are consistent within 1 <r, e.g. the temperature (metallic- 
ity) of the r < 3' region are 3.52 ± 0.17 (0.22 ± 0.07) and 
3.61 ±0.28 (0.22 ±0.11). 

In the image analysis for the surface brightness 
we carry out a vignetting correction for CDFS and 
CL0024+17, respectively. We derive the CDFS surface 
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Fig. 10. Left: Entropy (thick, black) and critical entropy floor (thin, black) using the global temperature (dashed) and using 
the observed temperature profile (solid, confidence intervals as dotted curves). Additional grey lines are the entropy floor (thin), 
S ~ 12AhjQ^ keV cm 2 , from Lloyd-Davies et al. (2000) and the predicted slope of 1.1 (thick) from the spherical accretion shock 
model (Kay 2004). The vertical line indicates the radius as defined in Fig. El Right: Scaled entropy of CL0024+17 vs. scaled 
radius using the radius of r2oo = 1.05/ifQ 1 Mpc (black, confidence intervals as dotted curves) and the Birmingham-CfA clusters 
in a temperature range of 2.9-4.6 keV (grey). 



brightness in the same detector coordinates as CL0024+17 
and scale it using the ctr ratio of CL0024+17 and CDFS 
in the 10-12 keV (12-14 keV) band for MOS (pn). 
After subtracting the scaled CDFS as a background from 
CL0024+17, we use a constant model to fit the residual 
soft X-ray background which is quite flat in the outer re- 
gion (5 < r < 10'). 

For comparison, we also used the blank sky back- 
ground from Lumb et al. (2002) and found that the pa- 
rameters change less than 1 a for both spectral and surface 
brightness analyses. 

2.3. Point source subtraction 

In a preliminary spectral analysis without subtracting the 
point sources we found contamination in both the soft and 
hard bands, especially for pn. 



Table 1. Parameters of the residual background model 
fitted in the 0.4-15 keV band. Col. (1): Instrument. Cols. 
(2-3): Index and normalization of the power law resid- 
ual background in which the normalization was scaled to 
1 arcmin 2 in units of 10~ 4 cts s _1 keV" 1 arcmin -2 . 



Instrument 


Index 


Normalization 


MOS1 


1.95 


1. 00 


MOS2 


2.21 


0.64 


pn 


1.95 


2.25 



Therefore we generated a list of the bright point 
sources using SAS task "edetect.chain" applied to five 
energy bands: 0.3-0.5 keV, 0.5-2 keV, 2-4.5 keV, 4.5- 
7.5 keV, 7.5-12 keV. The sources detected include fourteen 
of sixteen sources which were found in ROSAT HRI ob- 
servations (Soucail et al. 2000) except that S14 and S15 
are out of the FOV of XMM-Newton. Since we carried 
out analyses only in the r < 8' region, only the brightest 
fourteen point sources within this region (including nine 
ROSAT HRI sources, SI, S2, S3, S4, S5, S6, S8, S9 and 
S10 in Soucail et al. 2000) were subtracted from the source 
events using a radius of 40", which is comparable to the 
XMM-Newton Point Spread Function (PSF) cutoff radius 
of ~ 45", for each source that encloses nearly 90% of the 
flux of the point sources in both the spectroscopic and 
image analyses (cf. Table EJl. 

3. Results 

3.1. Image analysis 

The images in the 0.5-2 keV band for MOS and pn 
were corrected for vignetting using their weighted ex- 
posure maps and were then combined. We created an 
adaptively-smoothed image from the combined image with 
a maximum sigma (the width of the smoothing Gaussian) 
of 12" x 12" and a signal-to-noise of 5.0 (see Fig. P). 
Superposed XMM-Newton X-ray contours indicate an 
elongation in the northwest-southeast (NW-SE) direc- 
tion and substructure in the NW (R.A. = 00 ,l 26 m 25 ;i .8, 
decl. = 17°12'03".7, ~ 3.3' from the center). 




Fig. 11. Upper left: dark matter density and gas density, upper right: gravitational mass and gas mass, middle left: 

projected dark matter density (rwncate = 2.5/i^q 1 Mpc), middle right: projected gravitational mass (rtruncato = 2.5h^Q Mpc), 
lower left: projected dark matter density (rtruncato = 1.7h^ Mpc), lower right: projected gravitational mass (rtruncato = 
I.7/170 Mpc). The meaning of the lines is: (i) hydrostatic equilibrium (black) using the global temperature (thick, dashed) and 
using the observed temperature profile (thick, solid, confidence intervals as dotted curves) and its best fit (grey) of extended 
NFW model (thick) and NFW model (thin); (ii) HST lensing measurements (Kneib et al. 2003; dash-dotted). Additionally 
we present the gas density and gas mass (thin, black) with confidence intervals as dotted curves. We label the strong lensing 
projected gravitational mass of 1.59 ± 0.04 x W 14 h^M Q (Broadhurst et al. 2000) and 2.37 ± 0.36 x W 14 h^M Q (Tyson et al. 
1998) at the arc radii of 0.143/1^ Mpc and 0A53h^ Mpc and the Chandra results (Ota et al. 2004) using a /3/NFW model 
at the arc radius of 0.153/i^q 1 Mpc for comparison. The black vertical line indicates the radius as defined in Fig. El The g re Y 
vertical line indicates the arc radius of 0.143/if Mpc (e.g. Broadhurst et al. 2000). 
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We found no sign of the substructure described in 
Bonnet et al. (1994) which is ~ 6.9' from the cluster 
center (R.A. = 00 ,l 26 m 35 s .7, decl. = 17°09'35".8) to the 
northeast (NE). However, we observed the substructure, 
described in Kncib et al. (2003) and Czoske et al. (2001; 
2002), located ~ 3.3' NW from the cluster center covering 
1' with a temperature of 0.87 ± 0.13 keV fixing the metal- 
licity to Z = 0.3^0 and a bolometric X-ray luminosity 
of LY = 0.57 ± 0.06 x 10 43 /if 2 erg s" 1 , 2% of the total 
emission of the cluster. We also observed emission from a 
pair of background groups in the NE which were decribed 
in Czoske et al. (2002) at z = 0.495 based on a wide- 
field spectroscopic survey. The northern group centered 
at R.A. = 00 h 26 m 50M, decl. = 17°19'37".8 is more com- 
pact and X-ray luminous {L^ 1 = 0.24 x 10 44 /i^" 2 erg s _1 , 
r < 100"). The southern group centered at R.A. = 
00 ft 27 m 00 s .9, decl. = 17°14'51".7 is too faint for spec- 
tral analysis. If we assume a gas temperature in the range 
1-4 keV, then the bolometric luminosity of this group lies 
in the range 0.076-0.079 xlO 44 ^^ erg s" 1 (r < 100"). 

3.2. Chandra observation of the center 

CL0024+17 was also observed by Chandra (ID: 929). Ota 
et al. (2004) showed an adaptively smoothed image of the 
cluster with a rather symmetric appearance. 

A closer inspection reveals asymmetries, however. 
Fig.[2]shows an unsmoothed image (0.5-2.5 keV, ~2" pix- 
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Fig. 4. XMM-Newton HRM (~10" pixels) of CL0024+17. 
Additionally we show the position of the substructure (X- 
point) in kneib et al. (2003). 



els) and surface brightness profiles in annuli using the 0.5- 
2.5 keV image with azimuths 30°-120° (counter clockwise) 
from north to east and from south to west, respectively. 
There is a clear asymmetry between the northeast and 
southwest. To the northeast, we note a sharp decrease at 
a radius of about 34". The profile to the southwest dis- 
playing a more gradual decline. In Fig. |B| we show the 
azimuthally averaged Chandra temperature profile from a 
more detailed spectral study (kindly provided by Alexey 
Vikhlinin) fixing the metallicity to Z = 0.3-Zq. We note a 
corresponding temperature decrease approximately at the 
same radius as the surface brightness decline. This com- 
plex structure including a possible indication of a "shock 
front" might provide further support to the cluster merger 
scenario. 

3.3. Optical observations 

To compare the X-ray morphology with the optical im- 
age in the central region of CL0024+17, we retrieved 15 
HST WFPC2 observations including 11 from PI: Ellis (ID: 
8559) and 4 from PI: Turner (ID: 5453). They are block 
0, 11, 12, 13, 18, 19, 23, 24, 25, 36, 37 and 38 in Treu et 
al. (2003). All images were taken in the F814W filter (I 
band). We use the IRAF STSDAS and IMAGES packages 
to obtain the mosaic image from the calibrated science 
images. 

The X-ray contours of the XMM-Newton observations 
were plotted on the merged HST image (shown in Fig.|3l). 
We note the overlapping cluster centers and galaxy con- 
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Fig. 5. XMM-Newton spectra (top panel) extracted from the 
< 3' region fitted in the 1-10 keV band by an isothermal model 
and their residuals (bottom panel) for MOS (black) and pn 
(grey), respectively. 
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centrations in the NW (~ 3.3') of the X-ray and optical 
images. 

3.4. Hardness ratio maps 

To further study the substructure in CL0024+17 we pro- 
duced hardness ratio map (HRM) of the cluster, which 
primarily reflects the temperature distribution. The hard- 
ness ratio is the photon flux ratio of the hard band of 
2-7.5 keV to the soft band of 0.5-2 keV, HR = 1(2- 
7.5 keV)//(0.5-2 keV). An indication of an elongation in 
the NW-SE direction is visible in the HRM (see Fig. HJ) 
on a scale of r ~ 3.3'. The HRM shows a complicated 
structure with an average HR of 0.55 ±0.39 in the central 
r < 3' region. The substructure in the NW (~ 3.3') has an 
average of HR = 0.48±0.34 within 1'. We note two peaks: 
one is at (R.A. = 00 h 26 TO 24 s .8, deck = 17°11'13".2) 
that is marginally coincident with the NW substructure; 
the other is close to the center at R.A. = 00' I 26 m 30 s .9, 
deck = 17°09'50".3. 

3.5. Temperature and metallicity profiles 

We first derive global X-ray properties of CL0024+17 from 
the spectra extracted from the r < 3' region, covering radii 
up to a spherical overdensity of ~ 250, i.e. the ratio of the 
mean density of the dark halo with respect to the redshift- 
dependent critical density p cr jt(z). We measure an emis- 
sion weighted temperature of 3.52 ± 0.17 keV, metallicity 
of 0.22 ± 0.07 and bolometric luminosity of 2.9 ± 0.1 x 




0.1 1 

r (arcmin) 



Fig. 7. Flat fielded, point source subtracted and azimuthally 
averaged radial surface brightness profiles for CL0024+17 in 
the 0.5-2 keV band for pn (top) and MOS (bottom) and their 
best x 2 fits by a PSF convolved /3-model combining a con- 
stant residual soft X-ray background. Residuals scaled by the 
data uncertainties appear in the upper three panels for MOS1, 
MOS2 and pn, respectively, from top to bottom. 



10 44 /if 2 erg s- 1 (x 2 /dof = 248.5/228). We confirm the 
estimates of 3-4.5 keV and L<% 1 = 2.40 x 10 44 /if 2 erg s" 1 
in Bohringer et al. (2000) from the ROSAT observation 
and a estimate of 5.71 4 .' 2 keV in Soucail et al. (2000) 
from the ASCA observations. The Chandra observations 
yield a comparable temperature of 4.471q27 keV and a 
bolometric luminosity of 2.60 x 10 44 /i^ 2 erg s _1 (Ota et 
al. 2004). The source spectra for r < 3' and the fit- 
ted model are shown in Fig. [SJ A lower luminosity of 
1.26 x 10 44 /i7 2 erg s _1 is predicted from the L-T rela- 
tion in Arnaud & Evrard (1999) and the XMM-Newton 
temperature measurement. The high luminosity might in- 
dicate the complex cluster center which will be discussed 
in Sect.H 

For a more detailed spectral study we divide the clus- 
ter region into five annuli: 0-0.4', 0.4-0.8', 0.8-1.3', 1.3- 
2', and 2-3'. We investigate the dependence of the tem- 
perature and metallicity measurements on the low energy 
(low-E) cut-off (0.4 keV or 1 keV) as we did in our earlier 
study (Zhang et al. 2004a) to test the robustness of the 
spectral data fitting. We find that in the r < 1' region the 
derived temperatures increase by 15% if the low-E cut-off 
is increased from 0.4 keV to 1 keV. With the decreasing 
signal-noise-ratio (S/N) in the r > 1' region the derived 
temperatures increase by up to 40% if the low-E cut-off is 
increased from 0.4 keV to 1 keV. As in our earlier work 
(see also Pratt & Arnaud 2002) we attribute this change in 
the best fit temperature with energy range to the influence 
of a contaminating soft spectral component. To minimize 
its effect, we use the higher energy cut-off at 1 keV above 
which the results are stable. 

The radial metallicity distribution is very flat which 
is expected from simulations where merging efficiently 
flattens the metallicity profile (Kobayashi 2003). The 




0.5 r 



0.4 r i 

0.3 E i i i i = 

0.1 0.2 0.3 0.4 0.5 
r c (arcmin) 

Fig. 8. 90% confidence contours and the best fit values for 
the core radius r c and slope parameter j3 of the f3 model for 
the surface brightness of XMM-Newton for MOS1 (solid; tri- 
angle), MOS2 (dotted; square) and pn (dashed; diamond), re- 
spectively. Chandra (thin; Ota et al. 2004) and ROSAT (thick; 
Bohringer et al. 2000) measurements are also shown as crosses. 
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deprojection correction does not provide a significant 
change in the results. In the r < 1' (r > 1') re- 
gion the temperature measurements vary by a factor 
of 5% (50%), which is approximately within the er- 
ror bars. The temperature and metallicity profiles de- 
termined from the spectra in the 5 annuli (x 2 /dof = 
68.2/72, 77.0/78, 64.0/60, 50.2/60, 110.1/96 from the inner 
to the outer annulus) are shown in Fig. A variation of 
the temperature with radius is found here for the first time 
for CL0024+17. The inner three bins covering the radial 
range to a radius of 1.3' (416/ifQ 1 kpc) are consistent with 
an isothermal profile. This agrees with the Chandra re- 
sults by Ota et al. (2004) who found an isothermal tem- 
perature profile to an outer radius of 1.5' (480/1^ kpc). 
The form of the temperature profile is very similar to 
the profile we obtained for the sample of REFLEX-DXL 
(ROSAT-ESO Flux-Limited X-ray cluster survey, Distant 
X-ray Luminous) clusters (Zhang et al. 2004b). We scaled 
the temperature profile in Markevitch et al. (1998) to 
CL0024+17 using a radius of r 18Q = l.lQh^ Mpc (ob- 
tained in Sect. EH^) and an emission weighted global tem- 
perature of 3.52 keV. We find good overall agreement as 
found for the REFLEX-DXL clusters except that the flat 
part of the temperature profile extends to a larger radius 
than the average profile of Markevitch et al. (1998). We 
also scaled the temperature profile for Sersic 159—03 in 
Kaastra et al. (2001) to CL0024+17, and found a simi- 
larly rapid temperature drop as a function of radius as 
CL0024+17. 

We found an isothermal temperature of ~ 3.9 keV 
to a radius of 1.5' (480/iyQ 1 kpc) and a power law model 
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Fig. 9. Cooling time of the gas using the global temperature 
(dashed) and using the observed temperature profile (solid). 
The vertical line indicates r2oo = 1.7h^g Mpc, predicted from 
the lensing data in Kneib et al.(2003). 



with an index 7 = 0.98 outside this radius. We approxi- 
mate the temperature profile with two analytical models: 
(a) an isothermal model using the global temperature of 
3.52±0.17 keV; (b) a function with the shape of a Gaussian 
of k B T{r) = fceTg-e-lH^r 1 ) 2 ( see Table© which approx- 
imately fits the observed temperature profile. 

3.6. Surface brightness 

We obtained an azimuthally averaged surface brightness 
profile for CL0024+17 in the 0.5-2 keV band for each 
instrument, in which the central position is RA. = 
00 /l 26 m 35 s .7, deck = 17°09'35".8. This energy band 
is selected because it provides an almost temperature- 
independent X-ray emission coefficient over the expected 
temperature range. After the vignetting correction and 
double-step background subtraction described above, we 
find that a /3-model (e.g. Cavaliere & Fusco-Femiano 1976; 
Jones & Forman 1984) 

^(r) = S (l + 4r 3/3+1/2 C 1 ) 

' c 

convolved with the XMM-Newton PSF provides an ade- 
quate x 2 fit to the surface brightness profile (Fig. [7J|. A 
PSF at each position for each instrument is determined 
from an empirical calibration (Ghizzardi 2001) according 
to a mean energy ~ 1.25 keV of the 0.5-2 keV band and 
its off-axis radius. 

We compare the PSF convolved and unconvolved /3- 
modcl parameters in Table [21 and find that the core radii 
are overestimated and the slope parameters are under- 




Fig. 12. Gas mass fraction as a function of radius. The lines 
have an identical meaning to those in Fig. 1111 except that an 
additional horizontal line indicates the WMAP measurement 
fb = « b /fim = 0.166 (Spergel et al. 2003). 
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estimated in the latter case. The combined 90% confi- 
dence contours (see Fig. [SJ of the three instruments pro- 
vide a narrow range for r c and (3 (0.266 < r c < 0.295'; 
0.551 < (3 < 0.593) at the 90% confidence level. Bohringer 
et al. (2000) obtained r c = 0.11-0.28' and (3 = 0.425-0.550 
from the ROSAT HRI data corrected for PSF. Ota et al. 
(2004) obtained r c = 0.287 ± 0.014'and (3 = 0.55 ± 0.02 
from the Chandra observations. Our result confirms the 
results from previous X-ray observations with a higher 
accuracy because of the improvement of the image quality 
of XMM-Newton compared to previous observations. 

The surface brightness profile approximated by a (3- 
model can be analytically deprqjected to yield the emis- 
sion per volume element, £(r) = A(r)n^(r). With the given 
emissivity, A(r), from the applied plasma model, one can 
derive the electron density profile n c (r) = n e o(l + ^5-)~ 3/3 / 2 
with the parameters given in Table |3 combining the three 
detectors. We studied the influence of the central emission 
in the beta model analysis and found that the parameters 
like f3 and r c change within 2% if the central bins are 
masked. 

We obtain the pressure distribution P(r) as a function 
of radius directly from the temperature and electron den- 
sity measurements. Thus a /3-model is applied to fit the 
pressure with results shown in Table [S] 

Furthermore, we simulated a symmetric, background 
subtracted and flat fielded image using the parameters of 
the /3-model in Tabled The residual map, which contains 
the background, is obtained by extracting the simulated 
image from the adaptively-smoothed image. We confirm 
the existence of excess emission in the substructure in the 
NW, while the negative residual surface brightness in the 
cluster center shows that the elongation flattens the sur- 
face brightness compared to a symmetric structure. 



Table 2. Parameters for the brightest point sources de- 
tected in the r < 8' region. Col. (1): No. of the point 
source. Cols. (2-3): Central position. Col. (4): Flux deter- 
mined out to a radius of 20" in the 0.3-12 keV band. 



No. 


R.A. 


decl. 


Flux 








(10- 13 cgs) 


1 


00 26 12.8 


17 03 46.8 


1.317 


2 


00 26 30.9 


17 10 14.9 


0.933 


3 


00 27 00.0 


17 04 22.4 


0.890 


4 


00 26 31.0 


17 16 54.2 


0.760 


5 


00 26 20.1 


17 17 03.1 


0.677 


6 


00 27 05.3 


17 06 40.6 


0.671 


7 


00 26 44.2 


17 02 29.4 


0.581 


8 


00 26 26.1 


17 09 37.3 


0.558 


9 


00 27 07.5 


17 07 48.6 


0.542 


10 


00 26 17.2 


17 03 06.1 


0.400 


11 


00 26 45.8 


17 12 30.5 


0.391 


12 


00 26 17.9 


17 09 45.8 


0.259 


13 


00 27 03.7 


17 07 21.5 


0.224 


14 


00 26 42.8 


17 08 30.9 


0.060 



To further test if the cluster mass profile can be de- 
scribed by a NFW profile (Navarro et al. 1997; NFW), we 
fit the observed surface brightness profile by the model 
computed from the NFW model described dark matter 
halo (e.g. Makino et al. 1998). A model that fits the data 
well over a large range of radii shows an inner cusp in the 
surface brightness. The central surface brightness model 
at the given resolution is about a factor of two to three 
higher and thus inconsistent with our observations. Since 
we detect signatures of a cluster merger in the central 
region. We are not surprised to find deviations from the 
NFW model. Also Tyson et al. (1998) found a core in the 
mass profile in their lensing analysis. 

3. 7. Cooling time of the gas 

The cooling time is the total energy of the gas divided by 
the energy loss rate (e.g. Zhang & Wu 2003) 

°'v(t)(^) 4 (tf^)^< 2 » 

where n D and k^T are the electron number density and 
temperature, respectively, and g is the Gaunt factor (a 
function of k^T). The resulting cooling time as a function 
of radius is shown in Fig. |5| The central cooling time of 
4.5 Gyr is smaller than the age of the Universe and prob- 
ably smaller than the age of the cluster if we assume that 
the age of the cluster mass concentraion is longer than two 
crossing time because the merger frequency is expected to 
be a few Gyr (e.g. Schuecker et al. 2001). An even smaller 
central cooling time is implied by the Chandra data (Ota 
et al. 2004) due to the higher angular resolution of the 
central region. The data thus correspond, in the classical 
interpretation, to a small or moderate cooling flow (Fabian 
& Nulsen 1977). 

3.8. Gas entropy 

The main part of the observed entropy, defined as S = 
k^Tric 2 ^ 3 (e.g. Ponman et al. 1999), results from shock 
heating of the gas during cluster formation and it scales 
with the cluster critical temperature. An excess above 
this scaling law indicates the effect of an additional, 
non-gravitational heating source (e.g. Lloyd-Davies et al. 
2000). While a low central entropy indicates significant 
radiative cooling. We show the entropy profile derived 
for CL0024+17 from the XMM-Newton observations in 

Fig. ma 

The entropy of the r > O.lhjQ Mpc region derived 
from temperature measurements of CL0024+17 has a sim- 
ilar slope as that predicted from a spherical accretion 
shock model, S oc r (Kay 2004; thick grey line in 
Fig. I10p . The entropy under the assumption of isothermal- 
ity becomes steeper than the entropy derived from tem- 
perature measurements, especially in the outer region. 

The entropy of the gas in the center (r < O.lhjQ Mpc) 
lies below the entropy floor, S ~ 124/iyg 1 ^ 3 keV cm 2 , de- 
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Table 3. Parameters of the /3-model. Col. (1): Instrument. Col. (2): "Y" or "N" means the PSF convolved or 
unconvolved /3-model. Col. (3): Central surface brightness (0.5-2 keV) in cts s -1 arcmin~ 2 . Col. (4): Core radius in 
arcmin. Col. (5): Slope parameter. Col. (6): x 2 /dof . 





PSF 


So 


r c 


13 


x'/dof 


MOS1 


Y 


0.106 ±0.015 


0.27 ±0.03 


0.57 ±0.02 


99.7/63 




N 


0.047 ± 0.004 


0.34 ± 0.03 


0.56 ±0.01 


204.3/150 


MOS2 


Y 


0.094 ±0.010 


0.32 ± 0.03 


0.58 ±0.02 


79.2/65 




N 


0.045 ± 0.003 


0.36 ± 0.03 


0.54 ±0.01 


159.7/150 


pn 


Y 


0.375 ± 0.033 


0.25 ± 0.02 


0.56 ±0.01 


130.5/88 




N 


0.158 ±0.007 


0.33 ± 0.02 


0.54 ±0.01 


206.3/150 



rived for clusters by Lloyd-Davies et al. (2000). This indi- 
cates some effect from radiative cooling. 



We scale the entropy of CL0024±17 by 1 



fc B T(r) 
k B T 



where k B T — 2 keV is a constant related to the degree 
of preheating. In Fig. 1101 we compare it to the entropy of 
the Birmingham-CfA clusters in the temperature range of 
2.9-4.6 keV (Ponman et al. 2003) scaled also by l+ fc gjff . 
The entropy of CL0024±17 agrees with the self-similarly 
scaled entropy described by Ponman et al. (2003) within 
the observed dispersion of their cluster sample. 
Eq.J5J) can be rewritten as 



S r = lOOkeVcm 2 



2.869 x 10 10 yr 



f k B T 
VkeV 



.(3) 



If the cooling time t c is chosen approximately to be the 
age of the Universe, the above equation would allow us to 
determine the critical entropy floor S c of the gas in the 
cluster before the onset of significant cooling and possible 
mass deposition. In the center of CL0024±17 the observed 
entropy is slightly below this value (~ 154 keV cm 2 ) which 
implies a insignificant radiative cooling. 

3.9. Mass modeling 

We assume the intracluster gas to be in hydrostatic equi- 
librium within the gravitational potential dominated by 
dark matter. Neglecting the cosmological constant, A, at 
these high overdensities (a less than 1% effect, see Zhang 
et al. 2004a), we have 



1 d(n c k B T) GM(r) 



/im p r 



dr 



(4) 



where fi — 0.62 is the mean molecular weight per hydrogen 
atom. 

Analytic models of the gas density and pressure pro- 
files derived above can be easily combined with Eq.0} to 
obtain the mass profile (Fig. Illfl . The virial radius (not 
r 20o) is defined to be a radius with overdensity (the av- 
erage density with respect to p cr jt at the cosmic epoch 
z) of A c = 18tt 2 ± 82[ft M (z) - 1] - 39[0 M (z) - l] 2 for 
a flat universe. £Im(z) is the cosmic density parameter 
(e.g. Zhang & Wu 2003). Therefore the mass profiles in 
Fig. II II combined with the expression for the overdensity 
yield r V i r = 1.37/ifg 1 Mpc using the global temperature 



and r v ; r = l.26hjQ Mpc using the observed temperature 
distribution. 

We compare the XMM-Newton results with the re- 
cent measurements from Chandra and HST in Tableland 
Fig. The XMM-Newton measurements using the mea- 
sured temperature profile are slightly lower than Chandra 
values derived under the assumption of isothermality. As 
shown in Table 0] and Fig. ^2 the discrepancy between 
the X-ray measured total mass and strong lensing mass 
at the radii of 0.143/i 70 1 Mpc (Broadhurst et al. 2000) and 
0.153/ifo 1 Mpc (Tyson et al. 1998) remains large, about a 
factor of 4. 

Under the assumption of isothermality, the XMM- 
Newton measurements are more consistent with the ob- 
servationally determined M200-T relation, 

/ kT \ 1 - 60±0 - 04 
M 2nn = 4.11 ± 0.02 x 10 13 /if 1 ( ^ ) Mq , (5) 



VkeV J 



which is based on the conventional (3 model for the X- 
ray surface brightness profile and hydrostatic equilib- 
rium for 22 nearby clusters from Xu et al. (2001). For 
CL0024±17, Eq.© provides M 200 = 3.08 x 10 14 /i^Mq 
using the X-ray temperature. The masses obtained from 
the X-ray temperature measurement via the M-T rela- 
tions in Finoguenov et al. (2001) and in Bryan & Norman 
5.96 x 10 14 h^M e and M 200 = 
in which the former is more consis- 
tent with the HST optical lensing estimate. 

Navarro et al. (1997; NFW) describe a universal den- 
sity profile from numerical simulations in hierarchical clus- 
tering scenarios 

-1 / N -2 



ilTV 



Al 



2(H) 



(1998) 
3.76 x 10 14 h^M Q 



PDMM = ^critPcrit — 



1+^ 

r.. 



(6) 



where <5 C rit and r s are the characteristic density and scale 
of the halo, respectively, p C rit is the critical density of the 
Universe at the cosmic epoch z and p s = <5 C ritPcrit- #crit 
is related to the concentration parameter of a dark halo 
c = r vh ./r s by 

-1 



200 



ln(l + c) - 



1 



(7) 



We fit the X-ray determined mass profile with a NFW 
model, and found that the NFW model does not pro- 
vide a good fit to our data in the r < O.lh^ Mpc 
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Table 4. Comparison of the gravitational mass from the X-ray and optical lensing measurements. 



Observation XMM-Newton a Chandra " HST 

~7^ (h n l Mpc) Lil L05 O L7 

M 200 (lO^/i^Me) 2.3 ±0.1 2.0 ±0.3 4.6±g;^ 5.7^;J 

M proj (r = 0.143) (IO^/i^Mq) 0.61 ±0.01 0.48 ± 0.03 1.59 ± 0.04 



60 +0 14 

M P roj(r = 0.153) (10 14 ftf 1 M Q ) 0.64 ±0.01 0.51 ± 0.03 ' +0.13 [nfwi 2.37 ± 0.36 



0.53^111 (nfw) 
+0.7 



■1.1 



M(r = 1.4) (10 14 /if 1 M Q ) 3.0 ±0.2 2.4 ± 0.3 4.61 _ 
M(r = 1.7) (lO^/t^Mp) 3.6 ±0.2 2.8 ± 0.4 — 5^ 

a XMM-Newton results using the global temperature and using the observed temperature profile in which projected masses 
use the truncation radius of r 2 oo = 2.5h^ Mpc; 6 Ota et al. (2004); c Kneib et al. (2003); d Broadhurst et al. (2000); e Tyson 
et al. (1998). 



region (Fig. Ill|) . However the NFW model fit provides 
a virial radius of cr s = 1.28/iyQ 1 Mpc which is consis- 
tent with the virial radii of r vir = 1.37/iyQ 1 Mpc and 
r vir = 1.26/iyq 1 Mpc derived from the mass profiles us- 
ing the global temperature and using the observed tem- 
perature distribution, respectively, in Fig. 1111 combined 
with the expression for the overdensity. We derive a ra- 
dius of r2oo = 0-93hjQ Mpc from the empirical relation 
r 2Q0 = 0.813 (fcr B /keV) a5 (l±z)- 3 /2/ l7o i Mpc (Navarro et 
al. 1995), which is slightly lower than r2oo = l.ll/i™ Mpc 
using the global temperature and r2oo = 1.05/i^j Mpc 
using the observed temperature distribution. 

Furthermore, we apply an extended NFW model (e.g. 
Hernquist 1990; Zhao 1996; Moore et al. 1999), 



pDu(r) = <5critPcrit ( — 



r.. 



Q-3 



(8) 



to fit the data (Table which provides a negative index. 
This indicates that the dark matter density has a flat core 
for CL0024±17. 

We also present a comparison of the projected mass 
density and projected mass using the truncation radius of 

Table 5. Parameters of each profile model for the best % 2 
fit. 



Model 



parameter 



Gaussian 


k B T g (keV) 


4.2 ±0.4 




r g (ft^Mpc) 


0.08 ±0.04 




w 


0.69 ±0.04 




(ft^Mpc) 


0.093 ± 0.001 




n c0 (10~ 3 cm- 3 ) 


8.8 ±0.1 




/3 n " 


0.57 ±0.01 




r c p (fc^Mpc) 

Po(10- 2 keVcm~ 3 ) 

P P 


0.11 ±0.01 




3.59 ±0.04 




0.68 ±0.01 


NFW 


r s (ft^Mpc) 


0.37 ±0.03 




p s (10 14 M Q Mpc- 3 ) 


5.3 ±0.7 




c 


3.5 


extended NFW 


r s (ft^Mpc) 


0.072 ± 0.004 




Ps (10 16 M Q Mpc- 3 ) 


4.3 ± 1.7 




a 


-0.8 ±0.1 



either r2oo = 1-7/iyQ 1 Mpc, predicted from the lensing data 
in Kneib et al.(2003) or a larger radius of r = 2.5/itq 1 Mpc 
(Fig.lll)). The discrepancy is larger between the X-ray and 
strong lensing measurements than between the X-ray and 
weak lensing measurements. 

3.10. Gas mass fraction 

The gas mass fraction distribution according to the def- 
inition /gas(V) = M gas (r)/M(r) is shown in Fig. El We 

3/2 1 

found / gas = 0.19 ± 0.01/i 70 at r 2 oo = l.ll/i 70 Mpc us- 
ing the global temperature and 0.20±0.03fe 70 ' at r 2 oo = 
1.05/ifo 1 Mpc using the observed temperature profile, re- 
spectively. This result is comparable with the Chandra 

results / gas = 0.14to;o2 /l 7o 3/2 at r 2oo = lAOh^ Mpc (Ota 
et al. 2004), and agrees with the WMAP measured baryon 
fraction of the Universe /b = fib/^m = 0.166, where 
n h h 2 = 0.0224 and Q, m h 2 = 0.135 (Spergel et al. 2003). 
We obtain a gas mass fraction of 0.11 ± 0.01 using the 
global temperature at 7*2500 = 0.30 Mpc and 0.09±0.01 us- 
ing the observed temperature profile at r25oo = 0.33 Mpc, 
which agrees with the measurements of Allen et al. (2002) 
based on Chandra observations of seven clusters yielding 



0.105-0.138/i 



-3/2 

70 • 



However, this gas mass frac- 



tion at r2oo is slightly higher than the measurements of 
Sanderson et al. (2003) based on ASCA/GIS, ASCA/SIS 
and ROSAT/PSPC observations of 66 clusters yielding 

3/2 

/gas = 0.13 ± 0.01 h 70 , the measurements of Ettori et 
al. (2002) based on BeppoSAX observations of 22 nearby 
clusters, and the gas mass fraction for A1413 (Pratt & 
Arnaud 2002) at z = 0.143 based on XMM-Newton ob- 
servations yielding / gas ~ 0.12/i 



-3/2 

70 • 



The global parameters based on the XMM-Newton ob- 
servations for CL0024±17 are given in Table 

4. Discussion 

4.1. Metallicity 

We measure a global metallicity of 0.22 ± 0.07 Z Q which 
is typical for a cluster with an insignificant cooling flow 
(Fabian & Nulscn 1977). This metallicity value is in 
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good agreement with the averaged metallicity < Z >= 
0.2llo]o5 z o for 18 distant clusters with redshift 0.3 < 
z < 1.3 in Tozzi et al. (2003). However, Ota et al. (2004) 
obtain a much higher metallicity of 0.61-0.95 Z Q . We use 
the images in the 5-6 keV and 6-7 keV bands as indicators 
of the iron K line and continuum emission, respectively, 
and obtain an iron to continuum map. We extracted the 
spectra from a suggested high metallicity region in the 
iron-to-continuum map at R.A. = 00' I 26 m 32 s .7, decl. = 
17°09'46".3 covering a radius of 0.32'. No significant high 
metallicity was obtained by either setting a free temper- 
ature (5.4 ± 1.6 keV) or setting the temperature to the 
emission weighted temperature (4.09 keV) of the zone cov- 
ering radii of 0.4 < r < 0.8'. For the same radius of 
1.5' as used by Ota et al. (2004), we measure a metal- 
licity of 0.23 ± 0.07 and a temperature of 3.9 ± 0.2 keV 
(see Table 0) using the XMM-Newton data. Therefore the 
XMM-Newton data with higher significance provide no ex- 
planation for the high metallicity value obtained by Ota 
et al. (2004). The difference between XMM-Newton and 
Chandra metallicity measurements is an about 2 sigma 
effect, and may be due to the limited photon statistics in 
the Chandra data. 

4.2. Temperature gradient 

We derive an isothermal temperature of ~ 3.9 keV to a 
radius of 1.5' (480ft-f 1 kpc) and a power law model with 
an index of 7 = 0.98 outside this radius. The XMM- 
Newton temperature in the center is in good agreement 
with the Chandra measurement which extends to 1.5' in 
the Chandra observation of CL0024+17, covering the al- 
most isothermal center (r < 1.3') in the XMM-Newton 
observations. We provide a more extended temperature 
profile based on the XMM-Newton observations in which 
the temperature profile decreases in the 1.3-3' region. The 
strong temperature gradient at large radii makes the grav- 
itational mass lower by a factor of ~ 20-25% compared 
to the isothermal results. We found that the temperature 
gradient is not a sudden drop at a radius of 1.3' by mea- 
suring the temperature in different bin sizes (see Tabled). 
Within a radius of 3' for CL0024+17, we derived a poly- 
tropic index of 1.17 for the ICM in the outskirts which is 



lower than the adiabatic index of 5/3. This implies that 
the outskirts are still convectively stable. We consider the 
uncertainty in the mass estimate caused by the tempera- 
ture gradient by artificially increasing the slope parameter 
of the temperature distribution by a factor of ~20%, and 
find a decrease of the total mass by a factor of ~10%. 

Since the soft band (0.5-2 keV) is not so sensitive to 
the temperature map as the hard band (2-7.5 keV), the 
HRM can be used as an indicator of the temperature dis- 
tribution. Therefore the temperature map can be obtained 
from the HRM (e.g. Sanders & Fabian 2001; Sanders et 
al. 2003; Akimoto et al. 2003). The HR is also a reflection 
of absorption measure which affects the temperature (e.g. 
Akimoto et al. 2003). Thus we study the further proper- 
ties of the features in the HRM by performing a spectral 
analysis. For the statistic reason, we extracted the spec- 
tra from one of the peaks suggesting high temperature in 
the HRM at R.A. = 00 /, 26 m 30 <, .9, decl. = 17°09'50".3 
covering a radius of 0.33', and obtained a temperature 
of 5.8 ±1.0 keV with metallcity fixed at 0.3 Z© and 
5.7 ± 1.0 keV with the metallicity free {Z = 0.88±g;^| Z , 
i^ 01 = 0.43 x lO 44 /^ 2 erg s" 1 ). This supports the in- 
dication in the HRM of a high temperature region and 
not an artifact due to metallicity variations. However, the 
spectra can also be fitted by a power law with an index 
of 1.74 ± 0.14, which could then be attributed to some 
contamination from an AGN. This hot region is located 
exactly at the position of the substructure described in 
Czoske et al. (2001; 2002). 

4.3. Complex structure 

The Chandra image of the cluster core displays an indi- 
cation of complex structure. The projection of the counts 
as a function of radius using the high spatial rcsulotion 
image from Chandra shows an asymmetric structure and 
a surface brightness "edge" in the central region. The 
galaxy alignment in the HST image is perpendicular to 
the "edge" . Such an asymmetric central structure with an 
"edge" (r < 0.5') in the XMM-Newton HRM support the 
idea of a cluster merger. 



Table 6. Measured parameters for CL0024+17. Col. (1): 
Parameter. Cols. (2-3): Values obtained: (i) in the r < 3' 
region for line 2-4; (ii) using the global temperature and 
observed temperature distribution for line 5-7. 



Parameter Value 

kT B (keV) 3.52 ±0.17 

Z (Z ) 0.22 ±0.07 

L b £ l (10 44 h" 2 erg s" 1 ) 2.9 ± 0.1 

raoo (/ifo 1 M P C ) 1-11 L05 



M200 (lO 14 ftf o 1 M ) 2.3 ±0.1 2.0 ±0.3 
/ gas (fe TO 3/2 ) 0.19 ±0.01 0.20 ±0.03 



Table 7. Temperature and metallcity measurements in 
difference bin size of the radius. 



Region 


k B T (keV) 


Z (Z Q ) 


x'/dof 


r < 0.40' 


4.1 ±0.3 


0.16 ±0.10 


68.2/72 


r < 0.50' 


4.4 ±0.3 


0.10 ±0.09 


69.1/69 


r < 1.00' 


4.1 ±0.2 


0.25 ±0.08 


161.3/137 


r < 1.25' 


4.0 ±0.2 


0.23 ±0.07 


159.6/152 


r < 1.50' 


3.9 ±0.2 


0.23 ±0.07 


173.1/170 


r < 1.67' 


3.8 ±0.2 


0.24 ±0.07 


170.9/163 


r < 2.00' 


3.7 ±0.2 


0.23 ±0.07 


185.3/184 


r < 3.00' 


3.5 ±0.2 


0.22 ±0.07 


68.5/71 
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4.4. Comparison to lensing measurements 

The strongest disagreement between the X-ray mass pro- 
files and the strong lensing results is found in the clus- 
ter center. The difference between the X-ray determined 
mass (Table and the strong lensing mass of 1.59 ± 
0.04 x 10 14 /i^ M © at the arc radius of 0.143/ifg 1 Mpc 
(Broadhurst et al. 2000) and 2.37 ± 0.36 x lO^/i^M© 
at the arc radius of 0.153h^Q Mpc (Tyson et al. 1998) is 
up to a factor of ~ 4. 

This discrepancy is most probably due to the disturbed 
structure of the cluster center, highlighted by the elon- 
gated X-ray iso-surface-brightness contours in the central 
region. The substructure found in the gravitational lens- 
ing mass distribution map by Kneib et al. (2003) is also 
supported by the X-ray image. The velocity distribution 
of the galaxies can approximately be explained by a line- 
of-sight merger of two systems with a mass ratio of the 
order of 1:2 (Czoske et al. 2001; 2002). Therefore, the ap- 
plication of hydrostatic equilibrium to determine the mass 
in the central region of the cluster from the gas properties 
may not be valid. 

We observed an elongation in the NW-SE direction 
in the HRM. The negative residual surface brightness in 
the cluster center demonstrates that the elongation flat- 
tens the surface brightness compared to an azimuthally 
symmetric structure. The XMM-Newton image analysis 
shows substructure to the NW of the cluster center which 
is marked by a galaxy (see Fig. |3J in Kneib et al. (2003). 
This galaxy shows significantly extended X-ray emission 
compared to the XMM-Newton PSF and has a luminosity 
of L^ 1 = 0.57 ± 0.06 x 10 43 /if 2 erg s" 1 and a tempera- 
ture of 0.87 ± 0.13 keV within 1' at R.A. = 00' l 26 m 25 ;! .8, 
deck = 17°12'03".7. Since this high surface brightness 
region does not have a lower temperature but probably 
rather has a higher temperature than its surroundings, it 
does not mark a "cold front" but a high pressure region 
which coincident with a mass concentration seen in the 
lensing map. This implies that the second mass peak is 
included within the gaseous halo of the cluster and not 
only seen projected onto the cluster. 

4.5. Comparison to velocity dispersion measurement 

The M2oo _ cdm relation from simulations in Evrard & 
Gioia (2002) gives a comparable velocity dispersion of 
848^53 km s _1 from the HST mass and a relatively low 
value of the velocity dispersion of 598^^ km s" 1 using 
the XMM-Newton mass. Girardi & Mezzetti (2001) ob- 
tained a galaxy velocity dispersion of 911^®J 7 km s _1 
for CL0024+17. The velocity dispersion derived from the 
XMM-Newton mass is in good agreement with the veloc- 
ity dispersion measure in Czoske et al. (2002) who found 
velocity dispersions of 561^gg and 554lgg4 for the central 
component and foreground component which suggests a 
line-of-sight merger. If there is a line-of-sight merger, there 
are almost certainly additional galaxies in the line-of-sight 



which artificially increase the optical galaxy velocity dis- 
persion if only one component is assumed. 

4.6. X-ray properties 

The mass estimates from the M-T relation in Finoguenov 
et al. (2001) and in Bryan & Norman (1998) are M 20 o = 
5.96 x 10 14 hjjM Q and M 200 = 3.76 x 10 14 /i^Mg, 
which indicate a higher gravitational mass. We mea- 
sure a bolometric X-ray luminosity of 2.9 ± 0.1 x 
10 44 /i 

70 er g s 1 i which is higher than the luminosity 
of 1.26 x 10 44 h^Q erg s" 1 obtained from the L-T rela- 
tion in Arnaud & Evrard (1999) which might indicate a 
merger. The gas mass fraction at r 25 oo for CL0024+17 
agrees with the measurements of Allen et al. (2002) and 
Pratt & Arnaud (2002). The XMM-Newton measured gas 

mass fraction of 0.20±0.03ft-7 3/2 at r 20 o for CL0024+17 is 
higher than the values in Sanderson et al. (2003), Ettori 
et al. (2002) and Pratt & Arnaud (2002) for their clus- 
ter samples, which is possibly an indication for an under- 
estimated mass beyond r 2 soo as already noted from the 
M-T relation. Thus we suggest that the merger leads to 
distortions on global scales which reduces the mass for 
CL0024+17. However, on global scales the effect is smaller 
than in the cluster center. The relatively small distortion 
(less than a factor of 2) was found in several simulations 
(e.g. Evrard et al. 1996; Schindler 1996). 

4.7. Central cooling 

CL0024+17 has a bright central galaxy and a compact 
lensing core. The entropy in the center is a little below 
the entropy floor value suggested by Lloyd-Davies et al. 
(2000). However, the cluster almost shows little evidence 
for a cooling flow or a central temperature decrease which 
should be found in most cooling flow clusters. This sup- 
ports the scenario that the cluster center has been dis- 
turbed by a recent merger. 

5. Summary and Conclusions 

We performed a detailed imaging and spectroscopic study 
of the XMM-Newton observations of the lensing cluster 
CL0024+17, and obtained spatially resolved temperature, 
metallicity and density distributions. A temperature gra- 
dient was observed for the first time in CL0024+17 at 
large radii (1.3' < r < 3'). 

The image shows overlapping mass concentrations in 
the XMM-Newton X-ray data and HST optical data. The 
substructure which was observed by HST in Kneib et al. 
(2003) is confirmed in the X-ray data. The HRM shows 
an elongation in the NW-SE direction on a scale of 3.3'. 
Further spectroscopy confirms the temperature map sug- 
gesting high temperature region with k^T = 5.8±1.0 keV. 
Similar to Sersic 159—03, the observed temperature distri- 
bution shows a temperature gradient at 1.3' < r < 3'. The 
modeling of the temperature distribution yields a poly- 
tropic index of 1.17, lower than the adiabatic index of 
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7 = 5/3, which indicates a convectively stable state in the 
outskirts of CL0024+17. 

The cluster does not show a pronounced cooling flow in 
spite of a dominant central galaxy. Therefore we suggest 
that a merger in the cluster center (< O.l/i^ 1 Mpc) has 
disrupted any previous cooling flow. 

The mass estimate based on the precise measurements 
of the distributions of the temperature and gas density is 
lower by a factor of ~ 20-25% than the mass obtained un- 
der the assumption of isothermality in the outskirts. The 
NFW model does not fit the X-ray derived mass profile 
very well in the center (r < O.lh^Q Mpc). The extended 
NFW model fits the data but yields a negative index a. 
We attribute this to the disturbance of the core. 

The discrepancy remains between the X-ray gravita- 
tional mass and optical strong lensing mass at the radii of 
0.143/ifo 1 Mpc (Broadhurst et al. 2000) and 0.153h^ Mpc 
(Tyson et al. 1998) by a factor of up to 4. The XMM- 
Newton results are marginally consistent within 2a using 
the global temperature with the weak lensing measure- 
ment in Kneib et al. (2003). CL0024+17 is an example of 
a cluster for which weak lensing masses and X-ray masses 
are in acceptable agreement at large radii but for which 
the strong lensing masses disagree with the X-ray derived 
values at small radii (also e.g. A1689, Xue & Wu 2002; 
Dye et al. 2001; Clowe & Schneider). 

The XMM-Newton measured parameters for 
CL0024+17 deviate slightly from the empirical M— 
T and L-T relations. CL0024+17 has a bright central 
galaxy and a compact lensing core but does not show 
a significant cooling flow. Similar to A2218 (Girardi et 
al. 1997; Pratt et al. 2004), the apparent discrepancy 
between the X-ray and gravitational lensing determined 
mass of CL0024+17 is most probably due to a line-of- 
sight merger of two almost comparable subsystems. In 
addition, filamentary structures could also contribute 
to the projected mass detected by lensing, but since 
this mass does not lie within the cluster core, it would 
not be included in the X-ray mass measurement. The 
line-of-sight orientation makes it difficult to reveal the 
merger structure in the X-ray data, but on the other 
hand it enhances the probability of finding strong 
lensing features. Thus, the application of hydrostatic 
equilibrium assumption might break down in the central 
region. On global scales the agreement becomes better. 
The galaxy distribution, lensing mass maps, and X-ray 
data all provide a consistent description of the cluster 
morphology. 
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